Gravitational waves from galaxy encounters 
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We discuss the emission of gravitational radiation produced in encounters of dark matter galac- 
tic halos. To this aim we perform a number of numerical simulations of typical galaxy mergers, 
computing the associated gravitational radiation waveforms as well as the energy released in the 
processes. Our simulations yield dimensionless gravitational wave amplitudes of the order of 10~^^ 
and gravitational wave frequencies of the order of 10~^^ Ifz, when the galaxies are located at a 
distance of 10 Mpc. These values are of the same order as those arising in the gravitational radi- 
ation originated by strong variations of the gravitational field in the early Universe, and therefore, 
such gravitational waves cannot be directly observed by ground-based detectors. We discuss the 
feasibility of an indirect detection by means of the B-mode polarization of the Cosmic Microwave 
Background (CMB) induced by such waves. Our results show that the gravitational waves from 
encounters of dark matter galactic halos leave much too small an imprint on the CMB polarization 
to be actually observed with ongoing and future missions. 

PACS numbers: 04.30.Db, 98.80.Cq, 98.65.Fz, 95.35. -|-d 
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I. INTRODUCTION 



Galaxy encounters and interactions have well- 
established observational evidence. Supermassive black 
holes (SBHs) in the mass range of 10^ — 10^ Mq are firmly 
believed to exist at the cores of AGNs and quasars. Dur- 
ing such interactions gravitational radiation is going to 
be emitted. As a result, a stochastic background of grav- 
itational waves (GWs) from SBH encounters and coales- 
cence sweeps the cosmos, with wave frequencies of the 
order of 10^^ Hz or below. The gravitational waves pro- 
duced by the merger of binary SBHs, in terms of both the 
amplitude of the signal as well as the frequency range, are 
amenable to detection by the Laser Interferometer Space 
Antenna (LISA), becoming one of the most promising 
anticipated targets with an optimist signal-to-noise ratio 
of about 1000 

Binary black hole evolution and merger has received 
considerable numerical attention in recent years through 
N-body simulations (see, e.g. 0, Q and references 
therein). Examples include galaxy mergers in which the 
pre-merger galaxies contain power-law nuclear cusps and 
massive particles representing the SBHs. Previous work 
on the GWs emitted from such galaxy interactions have 
only considered the radiation produced during the coa- 
lescence of the central SBHs [j, Thus far, however, 
no attention has yet been paid to the gravitational radi- 
ation emitted from the bulk of matter of the interacting 
galaxies, i.e. neglecting the central SBHs, notwithstand- 
ing the fact that the mass of the galaxy which lies out- 
side the black hole is typically some orders of magnitude 
larger than the SBH mass itself. The underlying reason 
for the absence of investigations in such aspect of galaxy 
encounters may have to do with the fact that the size of 
the coalescing objects involved (the whole galaxies) is not 



compact enough to yield a direct observational measure 
in the frequency range of the upcoming interferometer 
LISA. 

On the other hand the dominant component account- 
ing for the total mass of the galaxy is the dark mat- 
ter part. Hence, an interesting question to ask oneself 
is whether a hypothetical detection of the gravitational 
radiation emitted from galaxy encounters, modelled as 
interactions of (coUisionless) dark matter halos lacking 
the SBH at the centre, would provide an unambiguous 
signature of the presence of such elusive matter. If that 
were the case the detection of those GWs would then 
provide a smoking-gun signature of yet another tradi- 
tionally evasive concept in gravitational physics, namely 
dark matter. 

Our aim in this work is to present and discuss numer- 
ical simulations of a reduced but representative sample 
of such galaxy encounters and to compute the associated 
gravitational waveforms and luminosities. Building on 
these results we conjecture about the detectability of such 
signals by electromagnetic means, namely through their 
imprint on the polarisation of the Cosmic Microwave 
Background (CMB hereafter). 

The non-linear evolution of cosmological scalar pertur- 
bations as sources of the CMB polarisation has been dis- 
cussed by [1] using perturbation theory. More recently Q 
have reported a related study which focus on the cosmo- 
logical stochastic GW background produced by idealised 
cold dark matter halos via power transfer from scalar 
and vector perturbations to tensor metric modes, during 
the strongly nonlinear stage of their evolution. Since the 
nonlinear evolution of such halos occurs on a cosmologi- 
cal timescale, the associated gravitational radiation may 
be relevant at frequencies comparable to those of the pri- 
mordial GW which affect the CMB photons and produce 
secondary CMB anisotropy and polarisation. While the 
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idea behind our work shares the views of 0, the GWs 
we analyse arise from galactic-hke halos of dark matter, 
which are substantially more asymmetric. Moreover, our 
study focus on the particular episode of galactic mergers, 
where the geometrical features of the process produce 
GWs with distinctive (burst-like) waveforms and larger 
amplitudes than those considered by fl] and therefore, 
their effects on the CMB polarisation are likely to have 
greater significance. 

Our study is further motivated by the prospects of po- 
tential experimental measurability. There are a num- 
ber of missions aimed to detect signatures of GWs (with 
wavelengths comparable to the size of the Universe) pro- 
duced by quantum fluctuations of spacctime during in- 
flation, by measuring the weak imprint they leave on 
the B-mode polarisation of the CMB. There have al- 
ready been claims Q suggesting that the polarisation 
induced on the CMB by relic gravitational waves at large 
scales can be already measured on the available data from 
WMAP, a situation which is likely to be improved by 
the Planck mission. According to [9| the B-mode po- 
larization spectrum of the CMB due to both, primor- 
dial (inflation) gravitational waves and also gravitational 
lensing, is well within the range of detectability of the 
Planck satellite. In addition, several experiments have 
been proposed for measuring the polarisation (E and B 
modes) of the CMB with higher sensitivity than WMAP 
(see e.g. [l^). Specially suited for the detection of B- 
mode polarisation induced by relic GWs are the upcom- 
ing projects CLOVER JT.] and BICEP [H. Further- 
more, among the missions included in NASA's Beyond 
Einstein program (see universe.nasa.gov) the Cosmic 
Inflation Probe also appears well suited to investigate the 
ultra-low frequency GWs from the early Universe. 

The simulations reported in this paper show that the 
frequencies of the gravitational waves produced by galaxy 
mergers (~ 10^^^ Hz) match the frequency sensitivity 
window of the aforementioned future missions (namely 
the Cosmic Inflation Probe). The amplitude of these 
GW signals is likely to be of comparable order than those 
coming from primordial waves originated during the early 
Universe. However, we estimate that the CMB polariza- 
tion they produce is weaker than that due to primor- 
dial waves. While the imprint of the gravitational radi- 
ation from galaxy mergers on the CMB polarization is 
inaccessible to detection, the gravitational lensing effect 
associated with gravity perturbations produced in our 
simulated galaxy mergers could marginally lie within the 
range of detectability of a number of planned missions 
designed to measure such polarization. This issue would 
be discussed elsewhere. 

The paper is organized as follows: In Section [H] we 
describe the galaxy model, the setup of the numerical 
simulations, and the numerical code used. In Section Hill 
the basic concepts and equations used to extract gravi- 
tational waveforms within the (Newtonian) slow-motion 
approximation are described. The results of our inves- 
tigation are presented in Section HVl Finally, Section fVl 



outlines the summary of this research. 



II. GALAXY MODEL AND NUMERICAL 
APPROACH 

Galaxy mergers are believed to be a formation mech- 
anism for present-day elliptical galaxies p^ . Different 
formation scenarios have been claimed to form high lu- 
minosity boxy ellipticals and low luminosity disky ellipti- 
cals. Boxy ellipticals could be formed through dry merg- 
ers of bulge dominated progenitors [13, [H, HI, while 
disky ellipticals could be formed through unequal mass 
merging of disks with a prominent bulge [18] or by in- 
cluding gas in the initial systems }19| . 

For the present study we adopt the coUisionless ap- 
proximation as a reasonable approach and simulate the 
merger of two elliptical galaxies and the merger of two 
disks with three different configurations (see Table [T|. 
The initial conditions for the models described here are 
taken from the samples studied by and We re- 
simulate three key experiments with a larger number of 
particles. We choose those experiments that produced 
the largest asymmetries during the merging stages for 
the remnants with realistic properties as compared to 
real-life ellipticals. 

To build and run the simulations we adopt non- 
dimensional units with G — 1 for Newton's constant of 
gravity. We use isotropic spherical Jaffe models 12011 as 
initial conditions for the elliptical galaxies (see [lj| for 
a detailed description of the algorithm). The projected 
surface density of such models exhibits a slope that de- 
creases with the radius roughly as R^^^, which makes 
it a suitable representation for elliptical galaxies. Jaffe 
model is characterised by the theoretical half mass radius 
rj, which we choose to be equal to 1. Because a cut-off 
radius is imposed at r = 10 x rj, the true half-mass radius 
ri/2 is 0.82. 

We follow [2l| to build the initial disk galaxy models. 
Our initial models have a bulge, a disk, and a halo. The 
initial parameters for these models are identical to models 
dbh in [l^ but with a larger number of particles. The 
model with a different mass is a scaled up version of the 
low mass system following a Tully-Fisher like relation 

For scaling our Jaffe models to elliptical galaxies the 
set of units is: [M] = Afj = 4 x 10" Mq; [L] = rj = 
10 kpc; [T] = 2.4 x lO"^ yr, and finally [v] = 414 km/s. 
A set of physical units that match the dbh models to the 
Milky Way are: [M] = 3.24 x 10" Mq; [L] = 14.0 kpc; 
[T] = 4.71 X 10^ yr, with [v] = 315 km/s. 

Table U reports the orbital parameters of the simu- 
lations. We first consider two equal-mass mergers be- 
tween elliptical galaxies with different impact parame- 
ters. Model Ml is a head-on collision on a parabolic 
orbit, while model M2 is a grazing encounter with a fair 
impact parameter, the orbital angular momentum is thus 
larger. Each elliptical galaxy is modelled with 4 x 10^ 
particles for model Ml and with 10^ particles for run 
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TABLE L Initial configurations for tlie models. (1) Model 
name, (2) type of merger, elliptical-elliptical (E+E) or disk- 
disk (S4-S), (3) Mass ratio, (4) initial orientation of the spin, 
(5) initial separation, (6) orbit ellipticity and (7) impact pa- 
rameter. 



Mod. 


Type 


Ml 
M2 


(6»2,02) 




e 


D 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Ml 


E-fE 


1:1 


(0,0);(0,0) 


40 


1 





M2 


E-fE 


1:1 


(0,0);(0,0) 


30 


1 


4.9 


M3 


S+S 


3:1 


(10,-10);(70,30) 


12.4 


0.7 


3.72 



M2. Model M3, on the other hand, is a merger between 
two disk galaxies with a mass ratio 3:1. Both systems 
are placed on penetrating prograde orbits. The orbits 
are elliptical (e = 0.7) to ensure merging and avoid large 
computational times. The two disks are placed with their 
spin slightly coupled with the orbital angular momentum. 
Each disk galaxy simulation consists of 2.75 x 10^ par- 
ticles for the halo, 1.75 x 10^ particles for the disk and 
9 X 10"* for the bulge. This makes a total of 10^ particles 
for the merger remnant. For all models of our sample, af- 
ter merger is completed the remnants are let to evolve in 
isolation for 10 half-mass radius crossing times to ensure 
that the inner parts are close to virial equilibrium. 

We have performed the simulations with the parallel 
version of GADGET [H using 16 CPU's on the lAC 
Beowulf Cluster. Each run takes of the order of 5 x 10^ 
seconds of CPU time. Independent softening is applied 
for each particle type. A tolerance parameter of 9 = 0.8, 
quadrupole terms, and a variable time step have been 
used. Energy conservation is satisfactory as variations 
are kept below the 0.5% level. 



III. GRAVITATIONAL RADIATION 

The study presented in this paper is focused in galaxy 
encounters taking place, in a flat Universe, at distances 
{D < 100 Mpc) much smaller than the horizon scale. 
The galactic systems we consider are far from being rel- 
ativistic in both senses, special relativity {v/c < 10~^), 
and general relativity {rg/R < 10~^, rg = 2GM/c^ being 
the Schwarzschild radius of the object). 

Under these conditions the gravitational radiation 
emitted from the encounters of our galaxy halos can be 
described by the so-called slow-motion formalism, where 
the spacetime metric can be linearised in the usual way 
(S/ij/ = fluv+h^v) and the transverse tracele ss ( TT) gauge 
can be used everywhere outside the system [2Jl ■ The spa- 

TT 

tial components of h^^^ are 

where lij are the components of the traceless inertial ten- 
sor. 



The collisionless particle contribution to lij is given by 
the following summation extended over all matter parti- 
cles. 




rip being the number of particles, TOp the mass of each 
particle, and 8ij Kronecker's delta. 

Outside the galaxies, the relative motion of two neigh- 
bouring particles A and B moving with the cosmologi- 
cal background (ideal detector) is fully described by the 

TT 

quantities /i^ . In the TT gauge there is a system of 
coordinates attached to A in which the coordinate varia- 
tions of the particle B are (t) - A^^, = iA^^ (t)]^ , 
where A' stands for the initial coordinates of the par- 

TT 

tide B and the quantities }l^ are calculated at point 
A. From this formula, it follows that oscillations in the 

TT 

h^j quantities lead to oscillations in the relative posi- 
tion of particles A and B with related frequencies and 
amplitudes. 

In the slow motion approximation, the gravitational 
luminosity, icWi is given by the well known formula 

A direct calculation of Law based on Eq.([3]) is prob- 
lematic [131 ^ ^ result of the noise introduced in 
the numerical computation of third-order time deriva- 
tives. However, as several authors have suggested (see, 
e.g. [lE[l^)i the second-order time derivatives involved 
in lij can be written in terms of the quantities x and i, 
which, in its turn, can be expressed in terms of related 
variables using the equations of motion of the system. 

It is well-known that, in the case of gravitational waves, 
the metric perturbation hij only depends on two inde- 
pendent linear polarisation states. Any electromagnetic 
radiation propagating in a spacetime perturbed from 
Minkowski by the presence of gravitational radiation, will 
suffer such polarisation itself. A particularly interesting 
case is represented by the CMB, and it is indeed pos- 
sible to take advantage of the polarisation of the CMB 
to design strategies to detect gravitational waves [1, [131 ■ 
The case of relic, ultra-low frequency GWs has been re- 
cently studied in great detail by [SJ . The order of magni- 
tude of the amplitude and frequency corresponding to the 
present-day primordial GW spectrum is ~ 10^^^ — 10~^ 
and - 10^2" - 10^1^, respectively (see Fig. 2 in Q). 

We will show in the following section that the source 
investigated in our work, galaxy encounters, leads to val- 
ues of the GW strain and frequencies which overlap to 
those arising in strong variations of the gravitational field 
in the early Universe. 
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FIG. 1: Dynamical evolution of the three merger simulations as seen from a point of view parallel to the orbital angular 
momentum. Left column: Ml, middle column: M2, right column: M3. Top row initial conditions, middle row right after the 
first pass through pericenter, bottom row at the end of the simulation. 



IV. RESULTS 

A. Dynamics of the encounters 

We begin the discussion of our results by analysing the 
morphological features which arise in our simulations of 
galaxy mergers, deferring for the next section the impli- 
cations of such dynamics in the emission of gravitational 
radiation. 

Figure [T] displays three snapshots of the three simu- 
lated merger encounters. Model Ml (Fig. [1] left column) 
is a head-on encounter between two equal-mass elliptical 
galaxies placed on an orbit with parabolic orbital en- 
ergy. The two galaxies approach each other and there 
is not much disturbance while they get closer prior to 
the first pass through the pericenter (which corresponds 
to the middle row of Fig. [T|), apart from a heating of 
the outer parts. A large fraction of the orbital energy is 



lost in the first encounter, a number of particles are ex- 
pelled from the system, becoming unbound and probably 
forming part of the unbound free floating intergalactic 
population [28]. These particles carry out a fraction of 
the orbital energy away leaving the cores of the merg- 
ing galaxies on a bound orbit. A rapid merging follows 
next, leading to the formation of a single object which is 
far from virial equilibrium right after the merging hap- 
pens. A reorganisation of the structure, orbital compo- 
sition of the stars and shape of the object takes place 
during some time, leaving behind a prolate non-rotating 
remnant when the simulation is stopped after the cen- 
tral parts of the remnant are close to virial equilibrium 
(bottom row of Fig. [1]). 

Model M2 (Fig. [T] middle column) is an equal-mass 
encounter between two elliptical galaxies, which follow 
a parabolic orbit with a large impact parameter. This 
simulation differs from the previous one in the type of 
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encounter. While for model Ml the encounter was head- 
on, model M2 corresponds to a grazing encounter. Again, 
prior to the first pass through the pericenter there is not 
much evolution in the two galaxies. After the first en- 
counter and due to the acting tidal forces, the galaxies 
feature extended plumes and broad bridges of particles 
between the two bodies. Like in model Ml, a fraction of 
the orbital energy is transfered to the particles as kinetic 
energy. A number of particles acquire positive binding 
energy and thus escape. After a second pass through 
the pericenter the galaxies finally merge. The orbital an- 
gular momentum of the initial galaxies is transferred as 
internal angular momentum to the stars of the merged 
system. As a consequence the final remnant presents a 
fair amount of rotation and due to the gentle merging 
orbit the final shape of the object is rather oblate. 

These two systems and merging processes give as a re- 
sult remnants that match the properties found in present- 
day, high luminosity elliptical galaxies. Therefore, this 
kind of dry-merger between elliptical galaxies has been 
proposed as a formation mechanism for present-day, high 
luminosity ellipticals [3, [IB] ■ 

Model M3 (Fig. [1] right column) corresponds to the 
collision between two disk-bulge-halo models with their 
spins partially aligned with the orbital angular momen- 
tum. Not much evolution of the disks is seen prior to the 
first encounter, although some distortion and heating is 
observed. However, after the first pass through the peri- 
center (middle row) prominent tails and a bridge between 
the two central bodies are created due to the spin-orbit 
coupling. The tidal tails carry away a number of stars 
that transport a large fraction of the angular momentum 
to the outer parts. After two more passes through the 
pericenter new tails are formed and finally the main bod- 
ies merge. Shells and ripples are still visible in the outer 
parts of the remnant after the inner parts have merged 
(bottom row) . The large bulges of the initial systems help 
to stabilise the disks during the merging stages, as dis- 
cussed by [III, so right after the last encounter and the 
merging a large fraction of the disk particles still have 
disk kinematics, leaving a clear signature in the merger 
remnant. Both the orbital angular momentum and the 
disks spin angular momentum are transferred to the in- 
ner parts of the remnant so the kinematics of this system 
is that of an oblate rotator and the isophotes of the rem- 
nant present disky deviations from perfect ellipses. This 
process has been proposed as a way to build present-day 
mid- and low-luminosity elliptical galaxies [3, HI] ■ 



B. Gravitational radiation: waveforms and 
energetics 

We turn now to describe the gravitational radiation 
generated in the galaxy encounters discussed in the pre- 
ceding section. Fig.[2]shows the dimensionless GW strain 
h (only the "-)-" polarisation is plotted) as a function of 
time for the three models. These waveforms have been 



computed using Eq. ([T]). This equation shows that the 
strain decays linearly with the distance. In order to com- 
pute the curves of Fig. [2] we have assumed a distance of 
D = 10 Mpc to the source, which is compatible with ob- 
servable galactic systems undergoing merger events (e.g. 
NGC 4424, The Antennae galaxies NGC 4038/9 or NGC 
520) located at few tens of Mpc. We have checked that 
the amplitude scales with the distance as expected, while 
the frequency remains unaffected. The mergers we are 
studying occur on characteristic timescales of a few 10^ 
years. Note that the radiation signal reaches its maxi- 
mum amplitude in such timescales. Therefore, such im- 
mense wavelengths lead to tiny GW frequencies of the 
order of - lO'^*^ Hz. 

From the very beginning of the simulations there is a 
non-negligible contribution to the GW signal associated 
with the internal dynamics of the halos. This is due to 
the fact that the halos are made of coUisionless particle 
distributions which are not in equilibrium and possess a 
nonzero quadrupole moment. All three waveforms are 
dominated by a maximum burst-like signal which is cor- 
related with the final plunge of the two galaxies, taking 
place at times ^ 2 x 10^ (top), ~ 1.8 x 10^ (middle), 
and ~ 3 X lO'^ (bottom) years, for models Ml, M2, and 
M3, respectively. Since the mergers extend considerably 
in time, even though the final snapshots in Fig. [1] show 
virialized systems (for their central parts), the waveforms 
keep memory of the violent events which leads to notice- 
able variations being still visible in the GW signal by the 
end of the simulations, particularly for model Ml. This 
effect, however, does not have any implication on the 
typical signal amplitudes {h ~ 10"^'^) and frequencies 
attained. 

The different dynamics of the mergers in models Ml 
and M2 is reflected in the waveform patterns shown in 
the upper top panels of Fig. [2l In both cases the plunge 
phase is clearly noticeable in the burst signal. Model Ml 
corresponds to a head-on collision, a process somehow 
more violent and shorter than that present in model M2, 
which describes a non-zero impact parameter encounter 
leading to a grazing collision. Despite the intrinsically 
different evolutionary tracks in models Ml and M2, the 
waveform gross features are similar. This could be ex- 
plained by the complexity of the scenario, and the way 
in which the waveforms are computed by adding the GW 
contribution from each dark matter particle, which may 
lead to a complex interference pattern. Therefore, the 
different dynamics in both models is hard to disentangle 
in the GW signal. 

On the other hand model M3 corresponds to a disk- 
like encounter. These system are inherently more as- 
pherical and, from the early phase of the evolution, the 
signal is stronger than in cases Ml and M2. As a direct 
consequence of the morphology of the disk-like galaxies 
involved in this encounter, tidal disruptions highly al- 
ter the disks for a more extended period of time than 
in the other two cases. Tidal streams are visible in the 
middle panel of Fig. [T]for model M3 (after the first pass 
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through the pericenter) , which translates in a modulated 
high-amplitude GW signal during most of the simulation. 
The signal peaks a,t t ^ 3.5 x 10^ years which, as in the 
other two cases, corresponds to the plunge event. Af- 
ter this, the signal decays keeping a high amplitude and 
variability due to the oblate geometry of the final object 
(see the spiral structure in the bottom panel of Fig. [T] for 
model M3). 

It is worth pointing out that since the values of h and 
the frequency of the GWs resulting from our simulated 
galaxy encounters are of the order of those expected from 
relic GWs (see e.g. Fig. 2 of their effect on the CMB 
polarisation is also likely to exist. 

We note that our results, in terms of both GW fre- 
quencies and amplitudes, differ from those of Q, being 
slightly larger for the case of galaxy encounters studied 
here. The explanation is to be found in the fact that Q 
study the formation and evolution of isolated, idealised, 
dark matter halos. No encounters are considered and 
the contribution to the GW signal comes only from the 
quiescence evolution of such halos. 

On the other hand, we have not included in our simu- 
lations the presence of SBHs in the numerical modelling. 
This has been done on purpose, motivated by the inter- 
est in searching for the contribution to the GW signal 
from the coUisionless component of the galaxies, leaving 
aside the contribution from the SBHs. This has no influ- 
ence for the signal we are studying since it is well-known 
that binary SBH mergers show chirp waveforms whose 
frequencies fall within the sensitivity range of LISA [ij 
and, hence, contribute in a extremely different band of 
the GW spectrum. 

Figure [3] shows the time evolution of the gravitational 
luminosity of the three models of our sample. The curves 
are computed using Eq. reducing the order of the 
time derivatives through the use of the equations of mo- 
tion. The first thing to mention is the correlation between 
the maxima of the luminosities and the GW amplitudes. 
As happens with the signal strain (see Fig. ^ the lumi- 
nosity in model M3 shows increased modulation than for 
models Ml and M2. We estimate that the average GW 
energy emitted during the whole duration of the simula- 
tions for all models is ^ 10'*^ erg. 

These results show that direct detection of GWs from 
galaxy mergers is entirely inaccessible to observation. 
In addition, we note that, unfortunately, gravitational 
waves from galaxy encounters cannot themselves polar- 
ize the CMB in a significant manner. An estimate based 
on ^d\, which is only valid for a source located in the 
wave zone (e.g. 100 Mpc), produces in our case a rota- 
tion angle of the polarization vector of the CMB photons 
^ 10^^'^ radians, a completely negligible effect. On the 
other hand, as mentioned in the introduction it has been 
recently claimed that the distribution of dark matter 
halos with masses between 5 x IO^Mq and IO^^Mq pro- 
duces a background of GWs whose effects on the CMB 
polarization deserve a detailed study. In particular, B- 
mode polarization could be detected in future missions 



tailored to study this polarization mode. This issue de- 
serves a detailed study which will be presented elsewhere. 



V. SUMMARY 

In this paper we have studied the gravitational waves 
produced by galaxy encounters, paying attention to the 
coUisionless component of the systems. To the best of 
our knowledge such an investigation has not been carried 
out before. For our study we have performed numerical 
simulations of a sample of representative galaxy encoun- 
ters, for which we have extracted, using the Newtonian 
quadrupole formula, the gravitational waveforms, ampli- 
tudes and frequencies, as well as the luminosities. The 
main result of our work has been to find out that GWs 
produced in galaxy encounters have typical values of the 
GW strain and frequencies which overlap to those arising 
in strong variations of the gravitational field in the early 
Universe. Our simulations confirm that the absence of 
the central SBH does not affect the GW spectra we have 
computed, since the signals from binary SBH mergers 
have very different features, in terms of frequency and 
amplitude, to the ones presented in this paper. 

Recently, a number of missions have been proposed 
to detect signatures of GWs (with wavelengths compa- 
rable to the size of the Universe) produced by quantum 
fiuctuations of spacetime during infiation, by measuring 
the weak imprint they leave on the polarization of the 
CMB. Given the fact that the main features of our re- 
sults could fit in the expected range of detectability of 
those missions, we have estimated the effects of a close 
galaxy encounter on the CMB polarization. Thus, we 
have investigated this possibility to find out whether this 
effect could add an unexpected source of GW "noise" 
which may hinder the original target of the experiments, 
namely relic gravitational waves from inflation. 

The most appealing motivation of this work could be 
the possible direct detection of dark matter in galaxy 
halos. By focusing on a particular galaxy encounter if 
the polarization of CMB photons were modified, this 
would imply that the source responsible of such effect 
would be the expected dark matter distribution. There- 
fore, dark matter should exist in concordance with the 
so-called dark matter paradigm. The detection of those 
elusive GWs (through their imprint on the CMB polar- 
ization) would provide a smoking-gun signature of the 
existence of the dark matter. However, the results of 
our estimates show that the gravitational-wave-induced 
imprints of galaxy encounters on the CMB are beyond 
observational capabilities. Nevertheless, there exist addi- 
tional aspects of the nonlinear dynamics of galaxy merg- 
ers which could have implications on potentially observ- 
able effects through CMB polarization, namely gravita- 
tional lensing [3l| or the so-called Stroskii effect (see [s^l 
and references therein). This issue will be discussed in a 
forthcoming paper. 

The results reported in this work are similar to those 
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investigated by ^3§\ and 0] regarding the nonlinear evolu- 
tion of cosmological structures as sources of gravitational 
radiation. Future work on this topic should account for 
the integrated effect of a realistic galaxy distribution. 
The hopes opened by the upcoming missions to mea- 
sure ultra-low frequency GWs make these investigations 
worth to pursue. 
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FIG. 2: Gravitational waveform for the three models: Ml 
(top), M2 (middle), and M3 (bottom). Only the "+" polar- 
isation mode is shown. Galaxies are located at f = 10 Mpc 
from the observer. 
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FIG. 3: Gravitational wave luminosity as a function of time 
for the three models: Ml (top), M2 (middle), and M3 (bot- 
tom). 



